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Abstract 

Let G be a simple undirected n-vertex graph with the characteristic polynomial of its 
Laplacian matrix L{G), det(A/ — L{G)) = ■ Laplacian-like energy of a 

graph is newly proposed graph invariant, defined as the sum of square roots of Laplacian 
eigenvalues. For bipartite graphs, the Laplacian-like energy coincides with the recently 
defined incidence energy IE{G) of a graph. In [D. Stevanovic, Laplacian-like energy of trees, 
MATCH Commun. Math. Comput. Chem. 61 (2009), 407-417.] the author introduced 
a partial ordering of graphs based on Laplacian coefficients. We point out that original 
proof was incorrect and illustrate the error on the example using Laplacian Estrada index. 
Furthermore, we found the inverse of Jacobian matrix with elements representing derivatives 
of symmetric polynomials of order n, and provide a corrected elementary proof of the fact: 
Let G and H be two n- vertex graphs; if for Laplacian coefficients holds Cfe(G) ^ Ck{H) for 
k ~ 1,2, . . . ,n — 1, then LEL{G) ^ LEL{H). In addition, we generalize this theorem and 
provide a necessary condition for functions that satisfy partial ordering based on Laplacian 
coefficients. 
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1 Introduction 



Let G = {V,E) be a simple undirected graph with n = \V\ vertices and m = \E\ edges. The 
Laplacian polynomial P{G, A) of G is the characteristic polynomial of its Laplacian matrix 
L{G) = D{G)-A{G), 

n 

P{G, A) = det(A/„ - L(G)) = ^(-l)fccfcA"-^ 

The Laplacian matrix L{G) has non-negative eigenvalues ^mi^ ■ ■ ■ ^ f^n-i ^ /^n = [Ij. 
Prom Viette's formulas, 

Cfc = 0-fe(^l,^2, • • • ,Mn-i) = ^ J|/^i (1) 

/C{l,2,...,n-1}, iGl 
\I\=k 

is a symmetric polynomial of order n — 1. Detailed introduction to graph Laplacians may be 
found in [2] and [3]. In particular, we have cq = 1, c„ = 0, ci = 2m, c„_i = nT{G), where t{G) 
denotes the number of spanning trees of G. If G is a tree, coefficient Cn-2 is equal to its Wiener 
index, which is a sum of distances between all pairs of vertices 

Cn-2{T) = W(T)= d(u,v). 

u,vGV 

The Wiener index is considered as one of the most used topological indices with high correlation 
with many physical and chemical properties of molecular compounds. For recent results and 
applications of Wiener index see [4] . 

Recently, Zhou and Gutman [5] proved that the extreme values of Laplacian coefficients 
among all n- vertex trees are attained on one side by the path P„, and on the other side by the 
star Sn- In other words, 

CkiSn) ^Ck{T) ^CkiPn), k = 0,1,2..., n, (2) 

holds for all trees T of order n. 

The Laplacian-like energy of graph G, LEL for short, is defined as follows: 

n-l 

LELiG) = Y,VJ^- 

k=l 

This concept was introduced in [6] where it was shown that it has similar features as molecular 
graph energy [7] (for recent results on graph energy see [8]). Various properties of Laplacian 
coefficients and Laplacian-like energy on trees and unicyclic graphs are investigated in [9j-|18j. 
It was demonstrated in [T2] that LEL can be used both in graph discriminating analysis and 
correlating studies for modeling a variety of physical and chemical properties and biological 
activities. 

The signless Laplacian matrix of the graph G is defined as Q{G) = D[G) + A{G). Matrix 
Q{G) also has all real and nonnegative eigenvalues n'l ^ ^'2 ^ ■ ■ ■ ^ fi'n ^ 0, [19]. Motivated by 
Nikiforov's idea [2D], Jooyandeh et al. [5T] introduced the concept of incidence energy IE{G) of 
a graph G, defining it as the sum of the singular values of the incidence matrix I{G). It turns 
out that 

n 

IE{G) = Y.^[. 

k=l 



In particular, if G is a bipartite graph, the spectra of Q{G) and L{G) coincide, and we have 
IE{G) = LEL{G). This relation provides a new interpretation of the Laplacian-like energy, 
and offers a new insight into its possible physical or chemical meaning. Many mathematical 
properties of this quantity were established in |22j and |23j . 

Stevanovic in ^lOj showed a connection between LEL and Laplacian coefficients. 

Theorem 1.1 Let G and H be two n-vertex graphs. If Ck{G) ^ Ck{H) for = 1, 2, . . . , re — 1 
then LEL{G) ^ LEL{H). Furthermore, if a strict inequality Ck{G) < Ck{H) holds for some 
1 ^ ^ re - 1, then LEL{G) < LEL{H). 



We will show that this fact remains true, but the original author's proof was not correct. 
Here, we will provide a corrected elementary proof. 

Let us move to a more general setting. Consider the open set in 

X = {(xi, X2, . . . , Xn) ■■ n > Xi > X2 > . . . > Xn-l > Xn > 0} . 

Let C denote the set of coefficients of polynomials having roots in X, 

C = {{ci,C2,...,Cn-l,Cn) ■■3{xi,X2,---,Xn) £ X 

P{x) = - C1X"~1 + C2x"-2 - . . . + (-l)"-lc„_iX + {-iTCn 
= {x- Xi){x - X2) ■ . . . ■ {x - Xn)] ■ 



Let F : X ^ C he the bijection defined by Viette's formulas ([T]). This function represents 
polynomial coefficients from C via the roots from X. Since every coordinate of the vector function 
F is a polynomial, we conclude that F is a continuous function. Symmetric polynomials are 
continuously differentiable functions and we have 



dxi 



E 

je/C{l,2,...,n-l}, 
\I\=k 



n ^i- 



(3) 



Let : C 



X be the inverse function of F. The Jacobian matrix of F is 

dci 



Jf 



dcj 
dxi 



i,j=l,n 



dci 
dxi 
dc2 
dxi 
dc3 
dxi 



dc„ 
dxi 



In [TU] the author has following identity 



<9ci 
8x2 
dc2 
8x2 
8c3 
8x2 



8c„ 
8x2 



8x„-i 

dc2 
8x„-i 

(9x„_i 



8x„- 



dci 
dx„ 
dc2 
dx„ 
dc3 
dx„ 



dcn 
dx„ 



dLEL 
dck 



EdLEL dfii 
1=1 «=i 



dLEL /dm 
dck/dfii 



(4) 



This is incorrect - we cannot use fact that the derivatives and |^ are reciprocal, since 
Cfc is a function of re variables. We can illustrate this error on the recently introduced molecular 
structure descriptor - Laplacian Estrada index [24j-|26]. defined as 



LEE{G) = ^( 



k=l 



Namely, the derivative ^gf^^ is strictly positive, since ^gf^ = e'^' > 0. Therefore, according 
to the equation this function is increasing in every coordinate Cfc. The following conclusion 
is obviously not true - from inequalities ([2D we have that LEL{Sn) < LEL{Pn) and on the other 
side 

n 



LEE{Pn) = Y,e 



2+2 cos ^ 



< e" + 1 + (n - 2)e = LEE{Sn), 



k=l 



since for n > 5 we have ne^ < e". 

This error can be fixed by finding the Jacobian matrix of the inverse F~^. Let 

Uj{x) = (x — Xl){x — X2) ■ . . . • {x — Xn) 

be the polynomial with n distinct real roots xi, X2, • • • , x„. The following theorem is the main 
contribution of this paper. 

Theorem 1.2 The Jacobian matrix of function F''^ equals 



J, 



(-1) 



n— 1 
LO'{Xi) 



- i,j=l,n 




(_-] \n-2 XI 

V ^> w'{xi) 

(_-] \n-2 X2 

V ^) iu'(x2) 
( i\n-2 x-j 

V ^' iv'ixa) 



-1 



1 



ui'(xi) 
(_-, \n-l 1 

1 



-11 



-2 Xri 



-1) 



n-1 1 



u}'{x„) 



The paper is organized as follows. After a few preliminary results involving Newton's iden- 
tities and complex analysis in Section 2, we prove Theorem 11.21 in Section 3. In Section 4 we 
provide a sufficient condition for establishing a partial ordering of graphs based on Laplacian 
coefficients and finally in Section 5 we give a corrected proof of Theorem 11.11 



2 Preliminary results 

Before presenting the proof of Theorem 11.21 ^6 need to prove some useful lemmas. 
Lemma 2.1 For every 2 ^ i ^ n and 1 ^ j ^ n holds: 

Proof. Notice that from the equation ^ , dci / dxj is a polynomial of degree i — 1 in variables 
xi,X2, . . . , Xj-i,Xj+i, . . . , Xn- On the right side, we have the sum of product of all {i — l)-tuples 
that does not contain Xj. We can derive this set if we consider all (i — l)-tuples and subtract 
those that contain xj - which gives the relation ([5]). □ 

Lemma 2.2 Let uj{x) be a polynomial of real coefficients of degree n having distinct real roots 
xi,X2, ■ ■ ■ ,Xn- For any nonnegative integer ^ k ^ n — 2 holds: 



Proof. Consider the rational function f{z) of complex variable z, defined as follows: 

By comparing degrees it holds |/(^;)| ^ -rfr for some constant C and large enough \z\. So, 
for the integral over the circumference Sr of radius R and center at the origin we have 

2ttC 



[ f{z)dz ^ [ \f{z)\dz^ [ ^dz 
Jsr Jsr Jsr 



R ' 



for large enough R. On the other hand, using the Cauchy Residue theorem we have: 



n n 

1=1 f{z)dz = 2m^Res{xi,f) = 27nJ2 



'Sr j^]^ 

Therefore, for large enough radius i?, we have inequality: 

C 



n u 

x7 



E 



R' 



and taking the limit as i? — t- oo we obtain that the sum equals zero. 



□ 



Lemma 2.3 Let uj{x) be a polynomial of real coefficients of degree n having distinct real roots 
xi,X2, . ■ . , Xn- It holds: 



n n-1 



1. 



Proof. Consider the function 



1 



uj{z) z zuj{z) 



Similarly, we get that \f{z)\ ^ -pj^. Using the additive property of integrals and 



dz r^"" iRe' 



z Jo Re^'t' 



the relation follows. 

For i ^ 1 denote by Si{xi,X2, ■ ■ ■ , Xn) the i-th power sum 



□ 



The Newton identities [28| (also known as Newton-Girard formulae) give relations between 
power sums and elementary symmetric polynomials. 

k 

k-Ck = ^(-l)'"icfc_i • Si. (6) 



i=l 



For A; > n, by definition is Ck = 0. We already proved that sq = si = . . . = 
Sn-i = 1. For A; = 2, . . . , n + 1, holds: 

{n + k - l)Cn+k-l = C„+fc-2Sl - Crt+fc-3S2 + • • • + (-l)"'~"^Cfc+iS„_2 

+ (-l)"-2c,S„_i + i-ir-'ck-lSn + ... + 

or equivalently 



Sn-2 = and 



CoSn+fc-l, 



Ck • Sn-1 — Cfc-1 • Sn + Cfc_2 • Sn+1 



■ bn+k-1 



0. (7) 



3 The proof of Theorem [Q] 



The theorem states that 

Jf ■ = Jp-^ ■ Jf = In = [^i,j\i,j=T^-: 

where I„ is identity matrix of order n. Therefore, we have to prove the following equalities for 
every 1 ^ < n, 

and 

Define the polynomials ooi{x), i = 1, 2, . . . , n as 

LOi{x) = {X - Xl){x - X2) ■ . . . ■ (X - Xi-l){x - Xj+l) ■ ... ■ {X - Xn). 

By the derivation of the product, we get 



n 

= y^^UJi{x) 
i=l 



The proof of ([8]) follows from the identity (case i = j): 

Uj'{Xi) = UJi{Xi) 

= {Xi - Xi){Xi - X2) ■ . . . ■ {Xi- Xi_i){xi - Xi+l) ■ . . . ■ {Xi - Xn) 

= ^ ■ xr' - ^ ■ xr' + p ■ xr' -■+ i-ir-'^ ■ + i-ir-'p. 

OXi OXi OXi OXi OXi 



For i ^ j, similarly we get: 



= u,{x,)=Y,{-l)'-' 



n~k 

k=l 



dxi 



Since the first row in Jp equals (1, 1, 1, ... , 1), from Lemma 12.21 and Lemma 12.31 follows that 
first row in the matrix product Jp ■ Jp-i is (1, 0, 0, . . . , 0). For the rows i = 2, 3, . . . , n we proceed 
by mathematical induction using the recurrent formula ([5]). 

n n—j " / f> \ 

EOCi _ / OCi-l \ Xf. 



^ dxk uj'{xk) ^ V dxk J uj'{xk) 



Ci-l 



E ^k _ \ ^ 



dci-i x\ 



\ ^'(xk) ^ dxk uj'{xk) 



For j > 1, we have that the sum equals —(—1)-' ""^(^i-ij-i = {—lySij. The case j = 1, we 
have to consider independently. 

" fir " f)r- , 

FIti i,i'('ri'\ j 1 n 1 ^/_^ 



^ dxk u}'{xk) ^ dxk uj'{xk) 

dxk uj'{xk)' 



EOCi-2 Xi 



Applying this substitutions i times and using the identity ([7]) we get that the first element 
in the i-th row equals zero. This finishes the inductive proof of Q. 



4 Partial ordering of graphs based on Laplacian coefficients 



First, we will give some general method for deriving the sign of special summation that we need. 
Let / be any smooth function and let 

P{x) = Cqx"""^ + Cix"~^ + . . . + Cn-2X + C„_i 

be the interpolating polynomial for / at the points xi,X2, ■ ■ ■ ,Xn of degree n — 1. Using the 
Lagrange interpolating theorem, we have: 

pf \_\^f( \ (X - Xl){x - X2) ■ {X- Xi-l)(x - Xj+l) ■ {X- Xn) 

i=i ^ ^' 

The sum 

/(x,) ^P{Xi) 



^ Uj'(Xi) ^ Uj'(Xi) 

i=l ^ ^' i=l ^ ^' 
equals to the leading coefficient of P. Using Taylor expansion of P{x), we have 

S = ±I^ = c, = ^-^^P^^^-^\x). 

Since f{x) — P{x) has n real roots, we conclude by Rolle's theorem that f^"'~'^\x) — P^"'~^\x) 
has at least one real root. Therefore, for some ^ in the interval containing all Xj's we have 

If the function /^"^^^(x) preserves the sign, we can uniquely determine the sign of sum S. 

5 The proof of Theorem 11.11 

Consider the open set in M"~-^ 

M = {(/ii,/i2, • ■ • : n > fli > H2 > . . . > Hn-l > 0} . 

Let C denote the set of coefficients of polynomials having roots in Ai, 
C = {(ci,C2, . . . ,c„_i) : (3(//i,/i2, • • • ,^n-i) e 

P(X) = - C1X"-1 + C2x"-2 + . . . + (-l)"-lc„_i 

= (x - /ii)(x - ^2) • • • (x - fJ-n^l)} . 

The Laplacian-like energy function LEL : — ?• M, defined by 

LEL{m,fi2, ■ ■ -,^71-1) = + \//^+ • • • + V/^n-i + 0, 

may then be represented as an implicit function of coefficients from C. By the chain rule, for 
arbitrary k, 1 ^ k ^ n, we have 

dLEL ^ dLEL dfn ^ 1 



^ dfii dck ^ 2^//iI dck ' 



Therefore, 



/77- ■ ^ ^ ■ ^/r,,.^ 9 



Let /(x) = x"'~'^~2 and according to the previous consideration we have to examine the sign 
of (n — 2)-th derivative of /. 

This means that the sign of /'■"^^^(x) is equals to (—1)'^""'^, and finahy 

dLEL 
-li — >0' 

OCk 

and the Laplacian-hke energy function is strictly increasing on C in each coordinate. 

So far we have dealt with the case of distinct eigenvalues only. The remaining step is to 
consider the closure of Ai 

M = {{pi,fi2, ■ ■ ■,fJ'n-l) : n ^ ^1 ^ ^2 ^ • • • ^ ^n-l ^ 0} . 

Its image under equations ([T]) is the set C of coefficients of polynomials having roots in Ai. 
The previous proof on Ai can not be applied since F is not injection - but, it can be done simply 
using the continuity. 

Suppose that we have two points xi and X2 in \ such that Cfc(xi) < Cfc(x2) for some 
k £ {1, 2, . . . , n — 1}, Cj(xi) = Cj(x2) for j k and LEL{xi) > LEL{x2)- Functions Ck and 
LEL are continuous, so we can find disjoint balls -Bi(xi,e) = {x | ||x — xi|| < e} and i?2(x2,e) = 
{x I ||x — X2II < e] such that Cfc(x) < Ck{y) and LEL{x) > LEL{y) for all x € Bi{xi,e) and 
y G B2{x2,e). The set \ has empty interior, so there exists points x' G Bi{xi,e) CiAi and 
y' G i?i(xi,e) n Al such that Cfc(x') < Ck{y') and LEL{x') > LEL{y'). Surfaces Cj{x) = Cj (for 
j ^ k and some constants Cj) can be arbitrarily close to the surfaces passing through points xi 
and X2 by continuity. Therefore, it follows that such surfaces have common points with both 
i?i(xi,e) and B2{x2,s) ~ so we can choose x' and y' such that Cj(x') = Cj{y') for every j ^ k. 

This is a contradiction with the first part of the proof, since x',y' G Ai. Finally, the 
Laplacian-like energy function LEL, which is strictly increasing on C in each coordinate, must 
be strictly increasing on C as well. 
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